INEQUALITIES FOR MODIFIED BESSEL FUNCTIONS AND 

THEIR INTEGRALS 



ROBERT GAUNT 



Abstract. Simple inequalities for some integrals involving the modified Bessel 
functions Ij,(x) and K v (x) are established. We also obtain a monotonicity re- 
sult for K v (x) and a new lower bound, that involves gamma functions, for 
K (x). 



1. Introduction and preliminary results 

In the developing Stein's method for Variance-Gamma distributions, Gaunt [5] 
required simple bounds, in terms of modified Bessel functions, for the integrals 



°JH u I u (t)dt and / e H u K v (t)dt, 



where x > 0, v > —1/2 and — 1 < f3 < 1. Closed form expressions for these 
integrals, in terms of modified Bessel functions and the modified Struve function 
L„(x), do in fact exist for the case (3 — 0. For z £ C and v E C, let S£ v {z) 
denote I v {z), e vm K v {z) or any linear combination of these functions, in which the 
coefficients are independent of v and z. From formula 10.43.2 of Olver et. al. [5] we 
have, for v ^ —1/2, 

(1.1) / z v Jf v (z)dz = ^T- x T(y + Ij2)z{^{z)-L v ^{z) - & v _ x {z)L v {z)). 



Whilst formula Q1.1JI holds for complex-valued z and v, throughout this paper 
we shall restrict our attention to the case of real- valued z and v. There are no 
closed form expressions in terms of modified Bessel and Struve functions in the 
literature for the integrals for the case j3 ^ 0. Moreover, even in the case j3 = 
the expression on the right-hand side of formula (11.11) is a complicated expression 
involving the modified Struve function L„(x). This provides the motivation for 
establishing simple bounds, in terms of modified Bessel functions, for the integrals 
defined in the first display. 

In this paper we establish, through the use of elementary properties of modified 
Bessel functions and straightforward calculations, simple bounds, that involve mod- 
ified Bessel functions, for the integrals given in the first display. Our bounds prove 
to be very useful when applied to calculations that arise in the study of Stein's 
method for Variance-Gamma distributions. We also obtain a monotonicity result 
and bound for the modified Bessel function of the second kind K v (x), as well as 
a simple but remarkably tight lower bound for Kq(x). These bounds are, again, 
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motivated by the need for such bounds in the study of Stein's method for Variance- 
Gamma distributions. Throughout this paper we make use of some elementary 
properties of modified Bessel functions and these are stated in the appendix. 



2. Inequalities for integrals involving modified Bessel functions 

Before presenting our first result concerning inequalities for integrals of modi- 
fied Bessel functions, we introduce some notation for the repeated integral of the 
function e l3x x L/ I v (x), which will be used in the following theorem. We define 
(2-1) 

Ifap.ofa) = e 0x x"I v (x), !{„,/}, n +i){x) = \ I(v,/3,n)(y)dy, 71 = 0,1,2,3,.... 



With this notation we have: 

Theorem 2.1. Let < 7 < 1> then the following inequalities hold for all x > 
(2.2) / t v I v {t)dt>x u I v+l {x), v>-l, 



(2.3) / t v I v {t)dt <x v l v (x) 1 v > 1/2, 



(2-4) /(„, ,n+l)(z) < l(v,0,n){x), V > 1/2, 

(2.5) J(„,- 7 ,n)(s) < )n e-^J K0 ,n)(a:), v > 1/2, n = 0, 1, 2, . . . , 

(2.6) / t v I u+n {t)dt < 2( f + n+ t % "I v+n+1 {x), i/ > -1/2, n > 0, 
' Iv + n + 1 



/ t v I u+n (t)dt 
Jo 

{" 2 
nf 
fc=i 

i r ™ 

< (1 _ 7) n (n 



2^ + 2fc 
2v + 2k 



x u I v + n {x), v > 0, n = 1,2, 3..., 



e-^a^Iv+nCa?), i/ > 1/2, n = 1,2,3,. 



2^ + fc 

Proof, (i) From the differentiation formula (|A.12I) we have that 

t v I u (t)dt= [ -t u+1 I L ,(t)dt>- f t u+1 I l ,{t)dt = x v I v+1 {x), 
Jo t x J 

since by (|A.2[) we have lim^o x v+1 I l/+ i(x) = for v > — 1. 
(ii) Using inequality (|A.7I) and then applying (|A.12[) we get 



t v I v (t)dt < / t v I v -i{t)dt = x u I v (x). 
Jo 

(iii) From inequality (|2.3[) . we have 

^,0,1) 0*0 < I{v,0fl){ x )- 

Integrating both sides of the above display n times with respect to x yields the 
desired inequality. 
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(iv) We prove the result by induction on n. The result is trivially true for n = 0. 
Suppose the result is true for n = k. From the inductive hypothesis we have 

(2.8) I {v _ 1>k+1) [x) = J I {v _ 1<k) {t)dt< , 1 _ nf s k J e lt h^o,k)( t ) dt - 
Integration by parts and inequality (|2.4[) gives 

e~ 7t I(„.o.k) (t)dt = e-"i x I {vfitk+1) (x) + 7 / e- 7 */ ( „. . fe+1) (t)dt 

Jo 

< e~ lx I [v ^ k+l) (x) + 7 / e-^I^o.k) (t)dt. 
Jo 

Rearranging we obtain 

e 7 I(u,o,k)(t)dt < YZT7f e 7X/ (":0,fe+i)(^), 
and substituting into (|2.8[) gives 

I{ v ,- 1>k+ i){x) < - +1 e~ lx I^ fiik+ i){x). 

Hence the result has been proved by induction. 

(v) From the differentiation formula (IA.12[) and identity (|A.10[) we get that 

j t (t»I u+n+1 (t)) = ±{t-^).t»+^I u+n+1 {t)) 

= t v I u+n {t) - (n + l)r- 1 / !/+ „+i(t) 

= - 2(p n + + n+lf Iu+n{t) + 2(Jl + n+lf Iu+n+2{t) 

_ 2v + n + \ n+1 
2(y + n + l) 2{y + n + 1) 

Integrating both sides over (0,x), applying the fundamental theorem of calculus 
and rearranging gives 

[ x tn v+n {t)dt = 2{ " + n+l K »i v+n+1 {x) - " + * 1 f ei v+n+2 {t)dt. 

Jo 2v + n+l 2v + n + 1 J 

The result now follows from the fact that l v (x) > for x > and by the positivity 
of the integral. 

(vi) From inequality (I2.6[) we have 

Jo 2v + l 

and 



I{v,o,2){x) = J I {ufl ,i){t)dt 
2v + l J 

2(^ + 1)2(^ + 2) VJ 

< — -— -x I v+2 {x). 

2v + l 2v + 2 + y ' 

Iterating gives the result. 

(vii) This follows from inequalities (|2.5[) and (|2.7p . □ 
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We now establish a simple lemma, which gives a monotonicity result for the ratio 



K j-^ . The lemma has an immediate corollary, which we will make use of in the 



proof of our next theorem. 

Lemma 2.2. Suppose x > 0, then the function K ^~j^ is strictly monotone in- 
creasing for v > 1/2, is constant for v = 1/2, and is strictly monotone decreasing 
for v < 1/2. 

Proof. To simplify the calculations, we let fi = u+1/2 and define h ll (x) — . 

It follows from the quotient rule and the differentiation formulas (|A.15|) and (|A.16j) 
that 

(2.9) h' M (x) = -l+^h li {x) + hl(x). 
Now, Theorem 2 of Segura [5] states that 

(2.10) fe»> /2 % = -M+V^ 2 + J^ fi>0iX>0 , 

and that for /i < the inequality is reversed and for = equality holds. Applying 
inequality (|2.10p to equation (|2.9p gives, for /i > 0, 



Similarly, we see that h' (x) = and that h'Ax) < for [i < 0. This completes the 
proof. □ 

Corollary 2.3. For v > 1/2 and a > 1 the equation K v (x) = aK v -i(x) has one 
root in the region x > 0. 

Proof. From the asymptotic formulas (|A.3[) and (|A.5[) . it follows that for v > 1/2, 

lim = 0, and lim = 1. 

xlO K v (x) x->oo K v [x) 

Since ^pjjff is strictly monotone decreasing on (0, oo), it follows that for a > 1 

the equation K v (x) — a.K v -\(x) (i.e. ^7p7~y^ — ^) has one root in the region 
x > 0. □ 

As an aside, we note that Lemma [2.3l allows us to easily establish an inequality for 
the Turanian A„(i) = K^(x) — K u _\(x)K u +i(x) (for more details on the Turanian 
A„(x) see Baricz pQ). 

Proposition 2.4. Suppose x > 0, then A u (x) < Aj,_ 1 (i) /or v > 1/2, A 1 / 2 (a;) = 
A_ 1 / 2 (x), and A„(x) > A 1/ _ 1 (.x) /or f < 1/2. 

Proof. By the quotient rule and differentiation formula (|A.14I) . we have 

d (K v -i(x)\ {K v {x) + K v _ 2 (x))K v (x) - {K v+1 {x) + K^x^K^x) 



dx \ K v {x) J 2K*(x) 

= gjj-iOjQ - K v _ 2 {x)K v {x) - {Kl{x) - K v ^{x)K v+1 {x)) 

2Kl{x) 

= A^a;)- A„(aQ 
2^) 
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Since, by Lemma 12.31 the function ^p7^ is strictly monotone increasing for v > 
1/2, is constant for v = 1/2, and is strictly monotone decreasing for v < 1/2, the 
result follows. □ 

With the aid of Corollary 12.31 and standard properties of the modified Bessel 
function K v (x), we can prove at the following theorem. 

Theorem 2.5. Let — 1 < < 1, then for all x > the following inequalities hold 
t v K v (t)dt < x v K v+1 {x), z/eK, 

(2.11) / t v K v (t)dt < x v K v {x), v<l/2, 



(2.12) / e^t v K u {t)dt < ——e 0x x v K v {x\ v < 1/2, 

J X ^ IP I 

(2.13) [°°1»K v (t)dt < ^ T{ L + l,2) x v K v {x), i/ > 1/2, 
Proof, (i) From the differentiation formula (|A.13|) we have that 

' pOO i i poo 

t u K w {t)dt= I -t v+1 K v {t)dt < - / t v ^ 1 K v {t)dt^x v K u+1 (x), 

J x t % J X 

since, by the asymptotic formula (IA.5I) . lim^^oo x" +1 K l/+ i(x) = 0. 

(ii) Using inequality (|A.8[) and then apply the differentiation formula (|A.13[) we 
have 



t u K„(t)dt< t v K v - X {t)dt = x v K v (x). 

J X 

(hi) Now suppose that v < 1/2 and f3 > 0. Using integration by parts and the 
differentiation formula (|A.13[) gives 

(2.14) / e fi H v K v {t)dt = --e fSx x v K v (x) + ~ I e 0t t" , K v - 1 {t)dt, 

J X P P J X 

Applying the inequality (|A.8|) and rearranging gives 

Inequality (|2.12[) for (3 > now follows on rearranging. 

The case (3 < is simple. Since e^* is a non increasing function of t when j3 < 
we have 

0t t v K v {t)dt < e 0x / t v K v (t)dt < e 0x x v K v {x) < —e 0x x u K v {x), 



where we used inequality (|2.11[) to obtain the the second inequality. Hence inequal- 
ity (|2.12[) has been proved. 

(iv) The case u = 1/2 is simple. Using (|A.1[) we may easily integrate t 1 ' 2 Ki/2{t). 

\^K 1J2 {t)dt = f°° J\e- f dt = J^e~ x = x^ 2 K 1/2 (x). 
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It therefore follows that inequality (|2.13|) holds for v = 1/2 because we have 

VET(i) _ 
r(i/2) ' 

where we used the facts that T(l) = 1 and T(l/2) = y/n. 

Now suppose v > 1/2. We begin by defining the function u(x) to be 



u(x) = Mx v K v {x) - / t v K v {t)dt, 

J X 

where _ 

^ + 1/2) 

r(iz) 

We now show that u(x) > for all x > 0, which will prove the result. We begin 
by noting that lim 2 ,_ ) .o+ u(x) = and lim^^oo u(x) — 0, which are verified by the 
following calculations, where we make use of the asymptotic formula (|A.3|) and the 
definite integral formula (|A.17[) . 

u(0) = lim ^-^±±^l x "K„(x)- (°° t v K v (t)dt 

= V^r( i . + l/2)2^ 1 - ^T{ V +\/2)2 v - 1 
= 0, 

and 

poo 

lim u(x) = lim Mx v KJx) - lim / t v K v (t)dt = 0, 

X— ¥00 X—^OO X— fOO J x 

where we used the asymptotic formula (|A.5[) to obtain the above equality. We may 
obtain an expression for the first derivative of u(x) by the use of the differentiation 
formula (|A.13I) as follows 

(2.15) u'(x) = x v [K v {x) - MK v -x{x)]. 

In the limit x — > + we have, by the asymptotic formula (|A.3|) . that 

V {2"- 1 r(i/)^ - M2\»-^Y{\v - 1|)^} , u^l, 
^ {2 v - 1 Y(v)^ 7 + M\ogx] , v=\. 

Since v > \v — 1| for v > 1/2 and lim a ._ > o+ x a \ogx = 0, where a > 0, we have 

u'(x) ~ 2 u ~ l T{v), asir^0 + , for v > 1/2. 

Therefore is initially an increasing function of x. In the limit x — > oo we have, 
by (EB, 



u'{x) 



u'{x) 



We therefore see that is an decreasing function of x for large, positive x. From 
the formula (|2.15|) we see that x* is a turning point of u(x) if and only if 

(2.16) K u[ x*) = ^^ K v _ l{ x*). 

From Corollary 12.31 it follows that equation (12.16[) has one root for v > 1/2 (for 
which ^ > 1). 
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Putting these results together, we see that u(x) is non-negative at the origin and 
initially increases until it reaches it maximum value at x* , it then decreases and 
tends to as x — ¥ oo. Therefore u(x) is non- negative for all x > when v > 1/2. 

(v) The proof for /3 < is easy and follows immediately from part (iv), since 1 < 
(i_ ) g2 2 )^+i/2 f° r v > 1/2. So we suppose 3 > 0. Again, because K 1 / 2 (x) = \f^ e ~ x ■, 
the case v = 1/2 is straightforward, so we also suppose v > 1/2. We make use of a 
similar argument to the one used in the proof of part (iv). We define the function 
v(x) to be 

/>oo 

v{x) = Ne Px x v K u {x) - / e pt t v K u (t)dt, 

J X 

where _ 

20^+1/2) 

(1- B 2 Y +1 l 2 T{u)' 

We now show that v(x) > for all x > 0, which will prove the result. We begin 
by noting that lim a ._ i .o+ v(x) > and lim. E _j. 00 v(x) — 0, which are verified by the 
following calculations, where we make use of the asymptotic formula (|A.3|) and the 
definite integral formula (|A.17|) . 



2"- 1 I»- / e^t v K v {t)dt 



2y^I> + l/2) 
(1 - 3 2 Y +1 / 2 T{v) 



(1 _ j83JH-l/2r( v ) 

0FI> + 1/2)2" 0Fl> + 1/2)2" 



0, 



and 



/>00 

lim u(z) = lim Ne Px x"K v {x) - lim / e pt t v K v (t)dt = 0, 

where we used the asymptotic formula (|A.5j) to obtain the above equality. We may 
obtain an expression for the first derivative of v(x) by the use of the differentiation 
formula (|A.13|) as follows 

(2.17) v'(x) = e? x x v [{l + Nf3)K u {x) - NK„-i(x)]. 

In the limit x —> + we have, by the asymptotic formula (|A.3|) . that 

' ' eP x x v {2"- 1 T(v){\ + NB)±j - N ■ 2^- 1 l- 1 r(|t/ - \\)-^} , v ^ 1, 
eP x x v {2"- 1 T(v)(l + NB)± + Nlogx} , v = 1. 



v'{x) 



As in part (iv), we see that v(x) is initially an increasing function of x. In the limit 
x — > oo we have 

t/(x) - (1 - JV(1 - B))^x"~ 1/2 e^- 1)x , for ^ > 1/2. 

Now, for v > 1/2 and < /3 < 1 we have, by (|A.5|) . 

x 2VttTV + 1/2) 1 1 1 

^ ^^)= i>) (i-^ ^ >M T L 
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Hence, v(x) is an decreasing function of x for large, positive x. From formula (|2.17|) 
we see that x* is a turning point of v(x) if and only if 

(2.19) (l+N0)K v {x*) =NK v -t(x*). 

Inequality (f2~18l) shows that N > 1 + N [3 for all v > 1/2 and < @ < 1. From 
Corollary |2.31 it follows that equation (I2.19[) has one root for positive x and therefore 
v{x) has one maximum which occurs at positive x. Putting these results together 
we see that v(x) is positive at the origin and initially increases until it reaches it 
maximum value at x* , it then decreases and tends to as x — > oo. Therefore v (x) 
is non- negative for all x > when v > 1/2, which completes the proof. □ 



Combining the inequalities of Theorems 12.11 and 12.51 and the indefinite integral 
formula (|1.1|) we may obtain lower and upper bounds for the quantity J£ v (x)\i v -\(x)- 
J£ v -\(x)Xu v {x) . Here is an example: 

Corollary 2.6. Suppose v > —1/2, then for all x > we have 

xV ~ H ^ <I I/ (x)L 1/ _ 1 (x)-I^ 1 (x)L 1/ (x) < ^ + 



V^-ir^ + 1/2) vy ' ' 1 vy ' ^2^-^(^ + 3/2) ' 

Proof. From the asymptotic formulas (IA.2I) and (|A.6[) for I v (x) and L(x), respec- 
tively, we have that 

\vav(x{I v {x)'L v ^i{x) — I u -i(x)~L u (x)) — 0, for v > —1/2. 

Therefore, applying the indefinite integral formula (11.11) gives, for v > —1/2, 

(2.20) f t v I v (t)dt = ^2"- 1 T{v + l/2)x{I v (x)h^ 1 {x) (x)U{x)). 
Jo 

From inequalities (|2.2p and (|2.6p of Theorem 12.11 we have 

x u I v+1 {x) < f t v I v (t)dt < 2 ^ + l K v I v+1 {x). 

JO Zl> + 1 

Substituting this inequality into (|2.20|) gives 

x v I u+x {x) < ^2 u - 1 T{v+l/2)x{I u {x)h u ^ 1 {x)- I u -!(x)L v (x)) 
2(f + l) ur i \ 

The desired inequality now follows from rearranging terms and an application of 
the standard formula xT(x) — T(x + 1). □ 

Remark 2.7. The lower and upper bounds for J ly (a;)Li / _ 1 (x) — 7 ly _ 1 (a;)L l/ (x) that 
are given in Corollary |2.6l are simple, but very tight for large v. 



3. Inequalities for the modified Bessel function of the second kind 

We now present some simple inequalities for the modified Bessel function of 
the second kind K v {x). The following theorem establishes an inequality for the 
modified Bessel function K v {x) that is useful in the study of Stein's method for 
Variance-Gamma distributions (see Gaunt [2]). 
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Theorem 3.1. Let fi > 1 and x > 0, then 

mi i x^k^x) 

1 • ' x 2 2»- 1 T{p) ' 

is a monotone decreasing function of x on (0,oo) and satisfies the following in- 
equality 

i t^~ 2 K (t) 1 
Proof. Applying the differentiation formula (IA.14I) gives 

d ( 1 x^- 2 K^{ x y 



(3.3) 



dx \x 2 2^" 1 r(/x) 

2 (/i - 2)^- 3 ^(a;) - \{K^ x (x) + K^ +1 {x))x^ 2 



x 



3 



2^'- 1 r(^) 



Using (|A.11[) we may simplify the numerator as follows 

1 

2 : 

= -a;if M _ 1 (s)-2if M (x). 
Hence, (|3.3j) simplifies to 

d / 1 x»- 2 K^{x)\ _ -2^T{^) + x^ +1 K^ l {x)+2x> 1 K^{x) 



(fi - 2)K ll {x) - -xiK^ix) + K M+1 (x)) 



= (jj, - 2)K fl (x) - -x 2K^ 1 (x) + -Z-K^x) 



dx \x 2 2^- 1 T( f i) J 2^- 1 T(^)x 3 
Thus, proving that (|3.1[) is monotone decreasing reduces to proving that, for x > 0, 
(3.4) x^ +1 K^-x{x) + 2x»K ll {x) < 2^T{ji). 

From (|A~13)l we get that 

^-(x^ 1 K^ 1 (x) + 2x^K^(x)) = ^-(x 2 ■x^ 1 K^ 1 (x)+2x^K^(x)) 

= 2x»K fi _ 1 (x) - x^ +1 K^ 2 (x) - 2x^K^ l {x) 
= -x» +1 K^ 2 (x) 
< 0. 

So x^ +1 K^i(x) + 2x fl K fl (x) is a monotone decreasing function of x and from the 
asymptotic formula (|A.3[) we see that its limit as x — > + is lim x _ > o+ (x^ 1 Kfj,-i{x) + 
2x»K ll {x)) = 2 ■ 2' 1 - 1 r(/i) = 2T(/i). Therefore (J33JI is proved, and so (^TJ) is 
monotone decreasing on (0, oo). It is therefore bounded above and below its values 
in the limits x — > oo and x — > 0. These are calculated using the asymptotic formulas 
(|A.5|) and (|A.4[) and are given below: 

lim 1 4-^3^1 = 0, 



z^oo \x 2 2^- 1 r( y u) 

1 x^- 2 K^{x)\ _ 2^- 3 r(^-i)_ 1 



lim — ^ — 



m+\x 2 2^- 1 T(p) J 2f- 1 T(fi) 4(^-1)' 

This completes the proof. □ 
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Finally, we establish a simple, but surprisingly tight, lower bound for the modi- 
fied Bessel function Kq(x). 

Theorem 3.2. Let x > 0, then 

Proof. Formula 10.32.8 of Olver et. al. 18! gives the following integral representation 
oiK (x): 



Setting t — 2u + 1 gives 



oo — 2xu 



K {x) = e~' / — du. 

Jo Vu + u 

For u > we have e 2u - 1 = > 2u + 2u2 > and 80 

/•oo ~2xu roo — (2s+l)« 

e x K (x) > V2 I - = = du = V2 I - 7 ===du, for x > 



Ve 211 - 1 7o Vl 3 



e 



Making the the change of variables y = e gives 

/•oo -(2x+l)« 1 /•! 

= -^=5(1/2,^ + 1/2) 

r(i/2)r(x + i/2) 

V2I> + 1) 
y/^T(x + 1/2) 
V2r(x + 1) ' 

where B(a, b) is the beta function, and we used the standard formula B(a, b) — 
^(l+b) ^° obtain the third equality. This completes the proof. □ 

Corollary 3.3. Let x > 0, then 

; < \ -e x K (x) < — . 

Proof. The upper bound follows because -ffofa) < Ki/ 2 {x) = \p^e~ x ■ The lower 
bound follows since r ^ x , +1 {? > , 1 which we now prove. Examining the proof 

r(x+l) y/x + l/2' K 6 -P 

of Theorem 13.21 we see that 

Tfx + 1/2) 2 f 00 e -( 2 ^+ 1 )« 

=<zu. 



r(x+l) 0Wo Vl-e" 2 "' 



Now, for u > we have 1 - e" 2 " = J2kLi (~l) fe+1 < 2u > and so 

r^+1/2) > 2 e-^)-^ = 2^ e _ (2x+i)v2dv = i 
r(a; + l) VWo V2w Jo Vd + 



y/x + 1/2 

as required □ 
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Remark 3.4. Luke [6] obtained the following bounds for Kq(x): 




Numerical experiments show that the bounds of Luke and our lower bound of Corol- 
lary [321 are remarkably accurate for all but very small x, for which the logarithmic 
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Appendix A. Elementary of properties modified Bessel functions 

Here we list standard properties of modified Bessel functions that are used 
throughout this paper. All these formulas can be found in Olver et al. [8], ex- 
cept for the inequalities and the integration formula (IA.17[) , which can be found in 
Gradshetyn and Ryzhik [3]. 

A.l. Basic properties. The modified Bessel functions I v {x) and K v {x) are both 
regular functions of x. They satisfy the following simple inequalities 




References 



I v {x) > for all x > 0, for v > — I, 



K u {x) > for all x > 0, for all v G R. 



A. 2. Spherical Bessel functions. 



(A.l) 
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A. 3. Asymptotic expansions. 

(A.2) I v {x)~ r(l / +1) (f) . x4.0,i/>-1, 

(A.3) ^(,)~/ 2|l/hlr ^l^ kl ' ^°.^ ' 

I — logo;, 2; 4. 0, v = 0, 

(A.4) jr„(a:) - 2 v-:l r(i/)ar" - 2 v-3 r(i/ - 1)^" +2 , 1 4. 0, 1/ > 1, 

(A.5) K v {x) ~ Jl^ e ~ X > x~^oo, 

(A.6) ^W~ ^ (|/ 2 +3/2 ) (f)' ,+1 ' ^0^>"V2. 

A.4. Inequalities. Let x > then following inequalities hold 
(A.7) I v {x) <I v -i (a), i/>1/2, 

(A.8) *r„(a!) < i^_x(aj), 1/ < 1/2, 

(A.9) ^(x) > ^-1(2:), ^ > 1/2. 



We have equality in (|A.9|) if and only if v = 1/2. The inequalities for K v (x) can 
be found in Ifantis and Siafarikas [4], whilst the inequality for l v {x) can be found 
in Jones [5] and Nasell [TJ. 

A.5. Identities. 



(A.10) /„+i(ar) = I v -x{x) L(x), 

x 

(A.ll) K v+1 {x) = K v - X (x) + —K v (x). 

x 

A.6. Differentiation. 

(A.12) ^-(x»I„(x))=x»I^ 1 (x), 
ax 

(A.13) jL( x » K „(x)) = -x v K v _ x (x\ 

ax 

(A. 14) = -^^+1^) + ^-lfc)). 

(A.15) ^{K v {x)) = -K v -i(x) - -K v (x), 

ax x 

(A.16) ^r{Kv{x)) = -K u+l {x) + -K v {x). 



A.7. Integration. 

(A.17) fj^K u (\t\)dt = ^^ffif , * > -1/2, -K /? < 1. 
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